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For /3-equilibrated nuclear matter we estimate the contribution to the bulk viscosity from purely 
leptonic processes, namely the conversion of electrons to and from muons. For oscillation frequencies 
in the kHz range, we find that this process provides the dominant contribution to the bulk viscosity 
when the temperature is well below the critical temperature for superconductivity or superfluidity 
of the nuclear matter. 

PACS numbers: 26.60.-c,21.65.Cd,67.10. Jn 



I. INTRODUCTION 



The bulk viscosity of nuclear matter plays an important role in the damping of oscillations in neutron stars. One 
well-known example is r-modes, which, if the interior of the star is a perfect (dissipationless) fluid, become unstable 
with respect to the emission of gravitational waves [IH3]. This emission acts as a brake on the rotation of the star. 
However, r-mode spindown will not occur if the r-mode is sufficiently strongly damped, for example by shear or bulk 
viscosity of the matter in the interior of the star. It is therefore important to calculate of the bulk viscosity of the 
various candidate phases in a neutron star. Several calculations exist in the literature, for nuclear [JHS] and hyperonic 
[I0l - [l2] as well as for unpaired quark matter [TSHIS] and various color-superconducting phases jT6l421j . 

In this paper we will study /3-equilibrated nuclear matter. We define the chemical potential for charged leptons 
to be /i; = — 0/e where (f) is the electrostatic potential and e is the positron charge. We will assume that the 
density is high enough that ^.i is greater than the mass of the muon, so the matter consists of neutrons, protons, 
electrons and muons. Such matter is expected to exist in the core of the star. In previous calculations of bulk 
viscosity of npe/i nuclear matter the focus has been on the contribution from interconversion of neutrons and protons 
via weak interactions. But nuclear matter at neutron-star densities is expected to show Cooper pairing of protons 
(superconductivity) or neutrons (superfluidity) either of which will suppress interconversion by a factor of 

order exp(— A/T), where A is the energy gap at the Fermi surface and T is the temperature. This opens up the 
possibility that, in superfluid or superconducting phases, the dominant contribution to the bulk viscosity might come 
from purely leptonic processes. The relevant process is conversion of electrons to muons (and vice versa) via either the 
direct Urea process or the modified Urea process. The direct Urea leptonic conversion process is forbidden by energy 
and momentum conservation: in converting an electron near its Fermi surface to a muon near its Fermi surface, the 
change in free energy is very small (of order T), so the emitted neutrinos carry momentum and energy of this order; but 

the change of momentum of the charged lepton is large, at least gmin — ^ii — ^ — , and the low-energy neutrino 

cannot carry this much momentum. However, the modified Urea process can occur; for example, two electrons with 
energy slightly above the Fermi energy can scatter to an electron and a muon with energies near the Fermi energy, or 
an electron and muon can scatter to two muons. The strongest interaction between leptons is electromagnetism, so 
this process proceeds via exchange of a photon, whose propagator should include the effects of screening by the nuclear 
medium. As the temperature decreases, the process will become suppressed as the Fermi distributions assume their 
zero-temperature step function profiles, but at finite temperature the modified Urea process will result in a non-zero 
contribution to the bulk viscosity. 

We calculate the leptonic bulk viscosity arising from the processes e + 1 ^ + £ + u + D, where t — e oi ji. All our 
calculations are in the "subthermal" regime where the density oscillation has a small amplitude, and the bulk viscosity 
is independent of that amplitude. We conclude that, if the protons and neutrons are both ungapped, i.e if there is 
neither superfluidity nor superconductivity, then the bulk viscosity from these purely leptonic processes is several 
orders of magnitude smaller than that from the nucleonic processes. However, once the temperature drops below the 
critical value for Cooper pairing of the protons or neutrons, the nucleonic bulk viscosity at frequencies > 10 Hz is 
strongly suppressed, and leptonic processes become the dominant source of bulk viscosity at those frequencies. 

In section |llj we lay out the process for calculating the bulk viscosity of a two-component leptonic system under 
application of a periodic volume and pressure perturbation. A crucial component of this calculation is the conversion 
rate between electrons and muons, which is discussed in HI In section IV we show the numerical results of our 
calculations and how they compare to the bulk viscosity resulting from modified Urea equilibration of the nucleon 
population. 
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II. BULK VISCOSITY OF LEPTONS 



First we write down a general expression for bulk viscosity in a two-species system, arising from interconversion of 
the two species. Then we specialize to the case of electrons and muons in nuclear matter. 



A. Bulk viscosity of a two-species system 



We assume that the system experiences a small-amplitude driving oscillation 

V{t) = V + Re{SV e"^^) 
p{t) = p + Re((5pe''"*) 



(1) 



where the volume amplitude 6V ^ ^ is real by convention, and the resultant pressure oscillation p[t) is complex. 
The average power dissipated per unit volume is 



dE 
~dt 

where t ~ 2tt/uj, so the bulk viscosity is [14] 



1 

7f 



2V^ dE 



Im((5p) V 
5V w 



(2) 



(3) 



We will determine Im((5p), which will be negative. We will assume that heat arising from dissipation is conducted 
away quickly, so the whole calculation is performed at constant temperature T. We assume that our system contains 
two particle species e and /i, and the state of the system is determined by the corresponding chemical potentials /ie 
and /i^. The total number of electrons and muons is conserved, and equilibrium is established via the conversion 
process e -f-)- /x. For simplicity of presentation and of the final expressions, it is better to work in terms of charged 
lepton number I and electron- muon asymmetry a, so pressure is a function of and where 



dp 
dm 

dp 

dfJ.a 



(4) 



From now on all partial derivatives with respect to fj,i will be assumed to be at constant fia, and vice versa. In 
beta-equilibrium, fXa is zero. The variations in the chemical potentials are expressed in terms of complex amplitudes 
Sfj.1, and Sfia, 



The pressure amplitude is then 



From ^ and ^ we find 



Sp — 



dp 
dm 



Sm + -7^ 
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(6) 



(7) 



To obtain the imaginary parts of the chemical potential amplitudes, we write down the rate of change of the corre- 
sponding conserved quantities. 



dni 


_ dni dpii 


dni dfia 


ni dV 


~di 


dp,i dt 


dfla dt 


~ ~V~dt ' 


dUa 


_ dua dpi 


^ dUa dfla 


_ ria dV 


dt 


dm dt 


dpa dt 


~V~dt ~ 



ntotal 
e— 



(8) 



3 



All the partial derivatives are evaluated at equilibrium, /x; = fli and fia = 0. The right hand term on the first line 
expresses the fact that charge is conserved, so when a volume is compressed, the density of charged leptons rises. On 
the second line, there is such a term from the compression of the existing population of particles, but there is also 
a rate of conversion F*'^'^ of electrons to muons, which reflects the fact that weak interactions will push the lepton 
densities towards their equilibrium value. For small deviations from equilibrium we expect F*'^'^ to be linear in fia, 
so it is convenient to write the rate in terms of an average width ja, which is defined in terms of the total rate by 
writing 



ptotal _ ^ (a\ 



We now substitute the assumed oscillations ([T]) and ^ in to and solve to obtain the amplitudes and 5^a 
in terms of the amplitude 5V and frequency lo of the driving oscillation. Inserting their imaginary parts in ([7| we 
obtain the bulk viscosity, which is conveniently expressed in terms of the susceptibilities 

dni 

dni dua 

Xla = 7^ = , (10) 

_ 

Ofla 

all evaluated at equilibrium, ni = fli, — 0. Note that Xai is the same as xia from Q. Defining 



7eff = 2" 7a 



nr^total 

XuXaa _ XU G'le^M 



c 



XllXaa - xla " XUXaa " xfa ^^a 
{Xll^a - Xlani)'^ 



(11) 



XiiixiiX 

we obtain the final result for the bulk viscosity in a two-species system, 



C = • (12) 



From ( 12 1 we can already see how the bulk viscosity of a two-species system depends on the frequency u> of the 
oscillation and the effective equilibration rate 7off. 

At fixed equilibration rate, the bulk viscosity decreases monotonically as the oscillation frequency rises; it is roughly 
constant for w < 7cff, and then drops off quickly as l/w^ for oj ^ 7cff. 

At fixed oscillation frequency the bulk viscosity is a non-monotonic function of the rate 7off. It is peaked at 
7off = with a value 

Cmax = iC/w . (13) 

For 7efi <C w or jeS ^ w the bulk viscosity tends to zero. Thus very fast and very slow processes are not an important 
source of bulk viscosity. As we will see below, for leptons in nuclear matter the equilibration rate is sensitive to 
temperature but the coefficient C is not, so we expect ({T) to be peaked at 7cff(r) — t^, where the oscillation 
frequency cj is of order kHz for typical oscillation modes of neutron stars. 



B. Leptons in nuclear matter 



In nuclear matter the leptonic chemical potential fj,i = = fJ-p, is much greater than the temperature and the 
electron mass, so we can evaluate the susceptibilities ( 10 ) at m,, — T = 0. Temperature dependence will come in only 
via the equilibration rate ja- Treating the electrons and muons as free fermions, we find 



7off 



C 



7a- 



(m; + kp? 
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9^ 



(14) 
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where the muon Fermi momentum is given by kp = fj,f — rn?^. Note that the bulk viscosity goes to zero as — 
(m^ — > me, reaUy). This is because if the muons and electrons have equal mass then under compression their relative 
densities do not change, and there is no need for any equilibrating process, so the pressure is always in phase with 
the volume and no dissipation occurs. 

Even without calculating the rate of lepton number equilibration, we can now estimate the amount of bulk viscosity 
that could possibly arise from leptons. If the equilibrating weak interaction at some temperature happened to 
have a rate that matched the typical oscillation frequency of the star, lo k, 2tt x 1000 Hz, and the lepton chemical 
potential had a relatively moderate value of about 120 MeV, we would obtain from (13), Cmax = 5.5 x 10^^ MeV"^ = 
7.5 X 10^^ gs^^cm"^. This is at the upper end of typical nuclear bulk viscosities which range up to lO'^^g s~^cm~^ 
[4]. This motivates us to proceed with the calculation of the rate of conversion of muons to and from electrons via 
the weak interaction. 

III. MUON-ELECTRON CONVERSION RATE 

The muon-electron conversion rate F^'^'^ consists of two partial rates, 

ptotal -ptotal 1 ptotal /'I '^'l 

e-^li ee— J-e/i e^-^fifj, \^) 

The partial rates are 

.total _ / d^Pi(Pp2(Ppzd?Pi(Pki(Pk2 



r^Uc. = j 64(2^)i4^i..2^3.;4l]il7, ^ iP^+P^-P^-P^-ki-k2)W,b^,a{piP2^P^P,k,k,) 

X [UiOi)h{L02) (l-/c(c^3)) (l-/<iK))-/cK)/d(c^2) (l-/a(^3)) (1-/0(^^4))] (16) 

where a,b,c,d are either e or fj,, Wab^cd is the spin-summed and averaged matrix element. The charged lepton of flavor 



j has energy ujj = yp'j +"^■^1 the neutrino of flavor j has energy ^lj = \kj\, and /h(wj) is the Fermi distribution 
function 



1 + exp 



(17) 



Using the previous definitions for /i; and fia, we have 

^ie^P^l+t^a, ^J'^^=^J■l-^J'a (18) 

and since is small, to first order in we have 



/e(c^l)/e(c.2) (l-/e(w3)) (1 " (^4)) " /e (^1 ) (^^2) (l-/e(a;3)) (l-/e(w4)) = F(c^l , ^^2 , C^3 , W4) ^ (19) 



and 



/^(C^l)/e(w2) (1-/m(w3)) (1-/p(c^4))-/^(^i)/m(c^2) (1-/^(^3)) (l-/e(^4)) = F(^l , 0^2 , C^3 , ^^4) ^ (20) 



2exp[(w3+W4-2/iz)/T][l + 2exp[(w2-AiO/7']+exp[(w2+^4-2M;)/7']] 
(l+cxp[(wi-/xO/r]) (l+exp[(L^2-A^0/7^]) {l+cMi^^- /T]) (l+exp[(c^4-Ai/)/T])' 

To determine the content of the matrix elements, we draw the Feynman diagrams for each possible way the reaction 
can occur. We can draw two different diagrams for each process, depending on the whether the weak conversion of 
the electron to muon occurs before the electromagnetic scattering, or in the reverse order (Fig. [ij Fig. [2|. However, 
because there are identical particles involved, and we are integrating over all initial and final momenta, we need to 
add two additional diagrams for each process. For the process e + e^/i + e + iz + P, we must add two diagrams where 
the labels on the initial state electron momenta are reversed; and for the process e + /i^ ^ + ^ + v + D, we must 
add two diagrams where the labels on the final state muon momenta are reversed. These diagrams get an additional 
negative sign for the interchange of fermions [2S]. For similar calculations, see [2^ 127) . 



FIG. 1: Feynman diagrams for the process e + e ^ e + n + + v^. There are an additional two diagrams which are obtained 
from these by exchanging pi o p2- 




FIG. 2: Feynman diagrams for the process e + /i— > fi + ji + Ufj, + Ve- There are an additional two diagrams which are obtained 
from these by exchanging o p4. 



Since we have four diagrams for each process, the spin summed-and-averaged matrix elements are 



Ml 

spins 



We^^^^ ^ l^^l + - Mg - (22) 

spins 



Here £'2, Es,E4 are the amplitudes corresponding to the diagrams of Fig. [T] and Mi, M2, M3, M4 are the ampli- 
tudes corresponding to the diagrams of Fig. [2] [55] : 

^1 = /of f ,M p^)re{pi)Pe{kih' (1-7^) /^ + ^+"'% 7pe(p2)M(p4)7A (1-7^) ^^(^2) 



E2 = J^^^/^^f iP^)re{pi)Ve{kl)l^ (1-7^) e(P2)/l(p4)7M ,,^.^_,ttr!:2 7A (1-7^) ^M(fc2) 

£^3 = ^. ^ —e 



H 

^0^^>3)7^eb2)Pe(fci)7^ (1-7^) (p;+!')r;2 7M^(Pi)A^(P4)7A (l-7^) -.(^2) 
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Ml = f^/ niPsWf^iPlPeiklh' (1-7^) /^;^^J"'% 7Me(P2)/i(p4)7A (l~7^) 



M2 = ' ,, M(P3)7V(Pi)^e(fci)7Ml-7^)^b2)Mb4)7M /" ^^^"^% 7A(l-7^)^^(fc2) 
M3 = "'f^ ,, A^(P4)7V(Pi)^e(fci)7^ (1-7^) /^^^,r''% 7Meb2)/i(P3)7A (1-7^) ^p(fc2) 
^^4 - "'f ^ ,, A^(P4)7Vbi)'^e(fci)7^ (1-7^) efe)A^(P3)7M /^"^J'"^ 7a (1-7^) ^.{k^) (24) 

where w = p2 — P3, and s = pi — p4. 

The only parameter in our calculation that depends on details of the baryonic matter in the neutron star is the 
plasma screening momentum Qs- In a full treatment one would have to use the appropriate in-medium propagator 
which is a complicated function of the photon momentum. 

In this paper we greatly simplify the calculation by assuming that the longitudinal and transverse photons have a 
common screening mass 

9? = 5a^if . (25) 

We argue in Appendix [X] that this leads to an estimate of the bulk viscosity that is correct to within an order 
of magnitude at reasonable densities for nuclear matter. As a further test we also performed calculations with no 
screening at all (g^ = 0) and found that the bulk viscosity shifted by no more than one order of magnitude. 



TO, 



To obtain the equilibration rates, we first multiply out the right hand sides of ( 22 ) and define partial matrix elements 

by 



y^lUe^ ^\Y. I^ll'' ^ie%e,, = ^ E (^1^2 + E\E{), W^^,^, = E ^^l^^ + ^3^^l)' ^te- (26) 

spins spins spins 

The traces resulting from the spin sums are easily evaluated with a computer algebra package; we used the FeynCalc 
package for Mathematica [29]. In the next few paragraphs, we will describe the steps used to analytically integrate 
10 of the 18 integrals, and list the expressions that we subsequently integrated numerically in Appendix [B} 

We make use of the fact that the neutrino energies are ~ T ^ /^e , /i^ by approximating the momentum and energy 
conserving delta functions as 

S^{Pl+P2-P3~P4~kl-k2) ~ 5(wi+W2-W3-t^4-f^l-f^2)^^(Pl+P2-P3-P4) • (27) 

We then note that ki and k2 occur exactly once in each term, dotted into one of the other 4- momenta pi . Writing 

kj ~ sin S^j cos r]j , sin sin rjj , cos ) (28) 

we can see that any dot product with another 4-momentum pi is 

Pi ■ kj = (w,— (pi)a; sin^j cos?7j-(pi)y sin^j sin?7j-(pi)2 cos^j) . (29) 
The integrals over the fci and /c2 angular variables then become trivial: 

^ ^ Pi ■ kj — J kjdkjd{coa ^j)drijPi ■ kj = Attloi J fl'^dflj (30) 

because all of the integrations over one of the angles S,j or rjj are zero. 

The energy-momentum conserving delta function allows us to use relations like Pi — p^ = Pa^ P2 to rewrite some 
of the denominators of the matrix elements. For example, in Wee^eii we can substitute variables so that p^ does 
not appear in the denominators of any of the terms; then we can integrate out the p3 3-momentum variables easily. 
Similarly, in Wefi-^^^ we can substitute variables so that p2 does not appear in the denominators and integrate out 
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the P2 3-nionientum variables. However, our matrix elements have many terms containing the four-momentum 
(P2), so it would be easier if we could integrate over d^pa (c?^p2)- This is accomplished by replacing 

d^Ps f d*p3 ( -\ f d'^ps 



, , ,-Si{pz),-M+K)^ / ^<5((p3)o-m0 (31) 

W3 J (P3)0 ^ ^ J [P3)0 

in Wee^eii and similarly for p2 in We^^^^ ■ In the last approximation we are using the fact that the Fermi distribution 
function is sharply peaked at low temperatures. Then we integrate over d'^pj, {d'^P2) using four of the delta functions. 

We can further approximate that the medium is isotropic, by taking one of the remaining momentum variables 
to be in a fixed direction (the z-axis for convenience). The electrons are relativistic, so uji — \pi\ and d^pi — 

ijjf duj id cos 9 id(j)i when particle i is an electron. The muons may not be relativistic, so uji — \/pf + rnf^ and d'^pi — 



I djf — m^dujid cos 9id(l>i when particle i is a muon. We then use the remaining delta function to integrate over the 

magnitude of this isotropic momentum variable. 

The remainder of the integrations are performed numerically. The only further approximation made was to again 
take advantage of the sharply peaked Fermi distribution function, and set uji = /x; everywhere inside the integral, 
except for inside the Fermi function itself. This allows a separation of the eight-dimensional integral into a four- 
dimensional energy integral and a four-dimensional integral over the angular variables. The integration variables are 
also changed to dimensionless variables by scaling them with respect to ni. 

The final expression for each term in the rate has the form 

„4^2 ,,12 , ,, , 

where £ is the species of the spectator lepton, and /f, and J^^ are dimensionless energy and angular integrals, 
respectively. These integrals are listed in Appendix [B| 

IV. NUMERICAL RESULTS AND CONCLUSIONS 

The remaining part of the rate calculation is performed numerically. The dimensionless energy integrals are nearly 
the same; a power-law fit of the results yields 

78.86 (;^)', - 78.62 g)' (33) 

In our approximation, the angular integrals only have dependence on /i;. We determined an analytical fit for the 
/i;-dependence of /^q and /^'q (accurate within 5%) over the range 120 MeV < /i; < 300 MeV by curve-fitting the 
numerical data with sixth-order polynomials: 

E^H « fi- 

1 .... . .,, ^ ^ 

Co = -1.7363 X 10^cl = 5.0189 x 1G^C2 = -4.7644 x 10'^, C3 = 1.3224 x 10^ 
C4 = 4.4203 X 10^C5 = -2.7199 x 10^, cq = 3.5119 x 10^ (34) 




Co 




1.2433 X 10^ ci = -3.6329 x 10^ C2 = 4.4365 x 10^ C3 = -2.8702 x 10^ 



C4 = 1.0354 X 10^C5 = -1.9728 X 10'\c6 = 1.5507 X lO'^ (35) 

Fig. |3] shows the fj,i dependence of the effective rate 7eff defined in ( 11 1. As approaches m^, the rate quickly drops 
to zero as the muon population disappears. The overall dependence is also illustrated in the sizable difference in 
order of magnitude of the rate for the three different temperatures. 

Fig. [4] shows the temperature dependence of the leptonic bulk viscosity ( as defined in (12), for an oscillation 
frequency w = 27r x IkHz. The three approximately straight lines on the log-log plot illustrate the power-law 
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FIG. 3: Dependence of the effective rate of electron/muon conversion 7eff (see (111) on tlie cliarged-lepton chemical potential 
fii at three different temperatures. As fii drops towards m^j, the muon population decreases and the conversion rate drops to 
zero. The temperature dependence is T^, hence j^fc is much larger at higher temperatures. 



dependence on T for three different values of . Also plotted are dotted curves showing the nucleonic bulk viscosity 
for two different values of the critical temperature Tc. These are obtained from Ref. [5] in a model where the neutrons 
are superfluid, pairing in the spin triplet state, the protons are superconducting, pairing in the spin singlet state, and 
they have a common critical temperature Tc — Tcp = Ten- Also, it is assumed that only modified Urea processes 
are available for damping of pulsations (although direct Urea processes would become possible at higher densities). 
Above the critical temperature for superfluidity/superfluidity, the bulk viscosity for 1 kHz oscillations due to leptons is 
several orders of magnitude less than the bulk viscosity due to nucleons. Below the critical temperature, the nucleonic 
bulk viscosity quickly decreases and at a low enough temperature, the leptonic contribution becomes dominant. Based 
on our calculations, this crossover temperature appears to be of order 0.01 to 0.1 MeV (10^ to 10^ K) for an oscillation 
frequency in the kHz range. Such a suppression of the nucleonic contribution can arise either from superfluidity of 
neutrons or from superconductivity of protons. It is therefore quite possible that for many cold neutron stars, the 
bulk viscosity of the superconducting or superfluid region comes mainly from leptonic processes. In regions that are 
neither superconducting nor superfluid (more strictly, where T > Tcp and T > Ten) the nucleonic bulk viscosity will 
likely dominate. 

The viscosity curves in Fig. [4] all slope upwards because the equilibration rate (T) is well below the oscillation 
frequency w, so we are in the slow-equilibration (high frequency) regime of (12), where 



7off(r) 



(36) 



This is true for both leptonic and nuclear viscosities. In this regime one can simply add the two bulk viscosities to 
get the total bulk viscosity (see, for example, appendix A of Ref. [17j). As the temperature rises, the equilibration 
rate and hence the bulk viscosity rise. When j^ff(T) comes close to ut, (36) becomes a poor approximation to (12): C 
reaches a maximum when "fcs(T) = lo. Those maxima, for both leptonic and nuclear bulk viscosities, are beyond the 
right hand limit of Fig. |4) for /x; = 200 MeV, the peak occurs at T « 40 MeV. 

We can now see how our results depend on the frequency of the oscillations. Decreasing lj moves each C,{T) curve 
to the left, shifting the viscosity curves in Fig. |4] upwards. The largest value we find for the leptonic effective rate 
(at r = 10 MeV, for /i; = 300 MeV) is 7eff ^ 2 rad/s, so for the leptonic bulk viscosity (36) is valid for oscillation 
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FIG. 4: (Color online) Dependence of the leptonic bulk viscosity on temperature for three different values of the lepton 
chemical potential, and an oscillation frequency of 1 kHz; for frequency dependence, see the discussion after ( 36 1 . We also show 
the nucleonic bulk viscosity [5] due to modified-Urca processes, for two values of the critical temperature. 



frequencies well above this value. For example, if we reduced the oscillation frequency from 1000 Hz to 100 Hz then 
all the viscosity curves in Fig. |4] would be shifted upwards by a factor of 100. Decreasing the frequency still further 
would bring us to the regime where, in the temperature range of interest, either the nuclear and leptonic rate was 
comparable to the oscillation frequency (so one or both bulk viscosity curves would show a resonant peak in our plot). 
Then one may not be able to simply add the bulk viscosities. At extremely low oscillation frequencies, both peaks 
would shift to very low temperatures, the bulk viscosity curves in our plot would all slope downwards, the nucleonic 
contribution dominates, and the bulk viscosities could again be added. 

It will be interesting to see whether the leptonic contribution that we have calculated here has any impact on 
oscillations of neutron stars. In the case of r-modes, shear viscosity becomes the dominant source of damping in the 
low temperature regime, so the leptonic contributions to the bulk viscosity at low temperature are not likely to be 
an important source of r-mode damping. Also the shear viscosity rj of superfluid nuclear matter is much larger than 
the leptonic bulk viscosity we have calculated: 77 ~ 10^^ gcm~-'^s~^ at T ~ 0.1 MeV, rising to 77 ^ lO'^^ gcm~^s~^ at 
T ^ 0.001 MeV (see Fig. 5 of Ref. [30]) so bulk viscosity would only dominate the damping of modes with very little 
shear flow. Radial pulsations [31(1321 would be an interesting example to investigate. We used a rough approximation 
(25) to treat the photon screening; we argued (Appendix |A| that this is valid to within about an order of magnitude, 
but if a more precise estimate of the bulk viscosity were required, one could improve on our treatment by replacing 
the approximation (25) with separate propagators for the transverse and longitudinal photons, incorporating their 
separate screening mechanisms |33j . It should be noted that our calculation is limited to the small-amplitude regime 
(Mo <C T). If the leptonic bulk viscosity is insufflcient to damp an unstable oscillation such as an r-mode then the 
amplitude will rise and it will be necessary to repeat our calculation in the large-amplitude ( "supra-thermal" ) regime 
[31] to see whether leptonic bulk viscosity can stop the growth of the mode once it reaches a large enough amplitude. 
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Appendix A: Photon screening 



In this appendix we discuss the adequacy of our approximation (25 1 for the internal photon propagator in the 
modified Urea process for leptons. The energy lo of the photon is ^ T because all the initial and final state particles 
have energies within T of their Fermi energies; however, the photon 3-momentum q must be large enough to move 
a lepton between the muon and electron Fermi surfaces, so g ^ qm\n where (jmin — PF,e — Pf,^,- Thus uj <^ q and 
we can write the photon self-energy in the static limit where it only depends on q. There are contributions to the 
longitudinal and transverse self-energies from protons, electrons, and muons. If the protons are superconducting, as 
they are at the temperatures of interest in this paper, then they provide an additional contribution to the transverse 
photon self-energy. The complete expressions are 



(Al) 

The Debye mass for a given species is (see, for example, [35]) 

Af|, = Ana^ikp (A2) 

where /i is the Fermi energy (defined relativistically, so = kp + rn^) and kp is the Fermi momentum. The screening 
functions and in the static limit are real, and are given by 

^2 , 



1 1 / r 

2 + 2^-(l-T 



) log 


q + 2 




q-2 



(A3) 

q + 2 



8V 4 q\ 4/ ^ q-2 
The full expressions for photon screening by a degenerate gas of charged fermions were first obtained by Lindhard 



[55] . Eq. (A3) was obtained from the version of Lindhard's expressions for the dielectric permittivities and Et 
given in Ref. [37], using the fact that ni(w,g) = (w^ — g^)(l — El) and UT{uj,q) = a;^(l — e^) (see Sec. (6.4) of 
Ref. [38 ). Note that and above are defined in the static limit, where w/g — > at fixed q. They are therefore 
different from the quantities x;(x) and Xt(a;) which are commonly given in the literature [33| I39j. and are calculated 

at uj = xq in the limit g — > 0. 

In our calculations of the leptonic flavor equilibration rate we use the rough approximation 11^= 11^ — qg (25) 



instead of the correct screening expressions given above. We now explain why this is a reasonable approximation. 
First we discuss the longitudinal photons. Th eir momentum varies from gmin up to kp^e + kp^fj_, but the momentum 



dependence of 11^ is very moderate: from (A3 1 we see that as q varies from to 2, varies from 1 to i. In order 
to judge whether, for the denominator of the longitudinal photon propagator, q^ -I- is a good approximation to 
q^ + Il]^{q), we use a naive free particle model for nuclear matter. In Table^we show the results. At each value of the 
baryon chemical potential fiB, the negative-charge chemical potential fii is determined by requiring overall electrical 
neutrality. This then fixes the Fermi momenta of the protons, electrons, and muons. In Table [ij we see that when 
q = Qmin (which is where there is greatest sensitivity to the exact form of the screening), the difference between q"^ + q^ 
and q^ + ni(g„iin) is a few percent at low density, and still less than a factor of 2 at very high densities. 

For the transverse photons, varies from at g = to | at g = oo, so the normal- fermion contribution to 
the transverse screening is more important at higher momenta. The other contribution to IIt comes from the 
superconducting protons, and it is more important at low momentum. At zero momentum we have Meissner screening, 
but as the momentum rises the effective screening mass drops slowly: this is seen in the calculation of Ref. [33] which 
finds that, for q 3> (where the correlation length ^ = pp^p/ {nipTc^p)), and assuming the static limit, 

<-^(g)^ ""y^^ . (A4) 

(This result follows from Ref. [33] eqn (49), taking w — >■ and using Q = tt^ as specified in the preceding paragraph.) 
In Table ^ we show numerical results for the naive free-nucleon model of nuclear matter. We assumed Tcp = 1 MeV 
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(see Ref. [23], fig. 10, and Ref. H^, fig. 2). At the lowest allowed photon momentum q = Qmin, which is where there 
is greatest sensitivity to the exact form of the screening, the difference between + and q^ + is a few percent 
at low density, but rises to a factor of 3 at density n/nsat = 12, and a factor of 10 at rUusat — 27. 

We conclude that our rough approximation of using a photon self-energy q^ — 5a^f ( 25 ) gives a reasonable estimate 
of the in-medium photon propagator. At low densities it is accurate to within 10%. At higher densities, up to 10 
times nuclear saturation density in the simple model of Table |lj our approximation underestimates the screening of 
longitudinal photons by a factor of about 2 and overestimates the screening of transverse photons by a factor of 
about 3. (At even higher densities, where a description in terms of nucleons is probably no longer appropriate, our 
approximation for transverse screening deviates further from the free nucleon model.) Since the rate involves the 
square of the photon propagator, we conclude that our approximate treatment of the photon propagator affects the 
rate by less than an order of magnitude at reasonable densities for nuclear matter. 



Ms 






2 

9min 


nL(<?min) 


nT(gmin) 


<?? = 5afif 


^min + nL((Jmin) 


9miii + nT(gmin) 


"Jmin + Qs 


1056 


3.164 


111.1 


5908 


1067 


55.45 


450.1 


6974 


5963 


6358 


1125 


6.76 


167.6 


1406 


2150 


30.13 


1025 


3557 


1436 


2431 


1200 


12.03 


224.4 


698.3 


3300 


65.06 


1838 


3999 


763.4 


2537 


1350 


26.93 


328.7 


304.3 


5783 


215 


3943 


6087 


519.3 


4247 



TABLE I: Screening parameters in MeV or MeV^ for a free-nucleon model of npefi nuclear matter; /is is the baryon number 
chemical potential, n/rigat is the baryon density relative to nuclear saturation density; /ij is the Fermi energy of the electrons 
and muons; Qmin is the lowest photon momentum that contributes to the modified Urea process; Hl and IIt are defined in ( Al I. 
The last three columns compare our approximate photon propagator at g = Qmin (final column) with the photon propagator 
using the full screening expressions given in Appendix \K\ 



Appendix B: Partial Rate Integrals 

The following abbreviations are used throughout this appendix: 



Cr2 = 

C24 = 

C'ls = 

Ci4 = 
F (Xq, X^, Xc, Xf/) 



Qs T 



1 — cos 6*2, Ci4 



x^ COS (J4 



1 — \J\ — x^^ (sin 6*2 sin Q4 (sin 02 sin 04 -I- cos 02 cos 04 ) -I- cos d-i cos O4 ) 
1 — (1 — x^„) (sin 01 sin 03 (sin 01 sin03-fcos0i cos 03)+ cos 0i cos 03) 

1 - (1 - X^) cos 01, ^34 = 1 - (1 - X^„)COS03 

2exp [(xc-Hxd-2) /t] [1-H2exp [(x^-l) A]-Hexp [(xfc-Hxd-2) A]] 

(l+exp[(xa-l)A]) (l+exp[(x6-l)A])' (l+exp[(xc-l)A]) (H-exp[(xrf- l)/t])' 



(Bl) 



= y dx2dx4dyidy2 yl yl F {x4 + yi+y2-X2 + l,X2,l,X4) 
I^^ = j dxidxzdyidy2 yl yl F [xi,l,xz,xi-X3-yi-y2 + l) 



(B2) 



jell 



dfl2dfl. 



4C12C14 + 2C12C24 — 2C14C24 — 4x^Ci2+X^C24 
^ (X2„-2C24-X2)2 



(B3) 



rel2 _ 

-'dn — 



dn2dfl4 



+ 



—2Ci2C24 + 2Ci4C|4-F 8(712(^14 -I-4C12C24 —4Ci2Ci4C24~4Ci4C24 
_ (X^-2C24-X2)2 

(4^422 -8C12+4C14+C12C24-C14C24) - xt 

(X2,-2C24-X2)2 



(B4) 



12 



^Ci2Ci4-4Ci2C24 + 6x^„Ci2 

\xl-2C2A-xj){xl-2C,4-x^,) 



(B5) 



7-el4 



+ 



2C12C14 — 4C12C14— 4Ci2C24 + 2Cl2C'i4C24 

(a;^-2C24-a;2)(ar^-2Ci4-a;2) 

, (-Ci^2 +6<gi2 - 2C12C14 +2Ci4 - 2C24) +x^/2 
(a;^-2C24-ar2)(x^-2Ci4-a;2) 



(B6) 



(ifi2<ii^. 



4Ci2Ci4+2Ci2C24-2Ci4C24-4x;'„Ci2+4x^Ci4+x^(724 - xt 
(x^-2C24-a;2)2 



(B7) 



re23 

-'da 



+ 



2C12CI4 — 4C12C14 — 4C12C24 + 2C12C14C24 

(x^-2C24-.T2)(xi-2Ci4-x2) 
(-^?2 + 6C12 - 2C12C24 - 2Ci4 + 2C24 ) +X^Ci2 /2 

(a.^-2C724-ar2)(:r^-2Ci4-a;^) 



(B8) 



Te24 



(if22<i^4 



+ 



2C'i2Cl4~4Ci2Ci4 + 2Cf2C'24~4Ci2C24 

(a;^-2C24-a;2)(a;^-2Ci4-ar2) 

^mHgl2 + 5gl2 + Cl2Cl4 + Ci2(724 + Cl4 + g24) 

(a;^-2C24-x2)(a;^-2Ci4-a;2) 



(B9) 



Te33 
-'dO 



2C12C14 +4C12C24 - 2C14C24 HC12 +C14) 
(x^-2Ci4-a;2)2 



(BIO) 



re34 

-'da 



+ 



-2Ci2C24+4Ci2Ci4 + 2C^4C24+8Ci2C24 ^4Ci2Ci4C24 — 4C14C24 
(X^-2C14-.X2)2 

rE^(2Ci^2-8gi2 + Cl2Cl4-Cl4g24+4C24) - 

(ar^-2Ci4-a;2)2 



(Bll) 



^do = / c«f^2rf^^. 



2Ci2Ci4+4Ci2C24-2Ci4C24+a;^HCi2+Ci4+4C724) - a;^ 



(a;^-2Ci4-ar2) 



2-\2 



(B12) 



-2(713^14 + 2(7i3(734 +4(714^34 (3Ci4 - 5(734) - a;^ 



(2x^-2Ci3-:r2)2 



(B13) 



jij.11 _ 
-'da — 



-{1-xlf^ J dQidns 



-2Ci3C^4 + 4C^4 — 2C13CI4 — 4CI4 + 8C13C34 + 8C14C34 

(2a;^-2(7i3-ar2)2 



+ 



+ 



^to(4Ci4 + 4C34 — 8C13 + 2C13C14 — 8C14 — 2C13C34 — 4C14C34 — 8C34) 

{2xl-2C,s-xir 
xf^ (3(7i3 - 2Ci4 + 2C34 + 10) - 3a:^ 



(B14) 



13 



T/il3 



-(1 



X. 



2 \3/2 



dQ.idQ.z 



2C13CI4 + 2C14CI4 — 4C13C34 — 4C14C34 



+ 



+ 



(2x^-2Ci3-x2)(2x^-2Ci4-x2) 

a^m Hg|4 + 5Cl3 + 5Ci4 + C13C34 + ^14^34 + 5(734) 

(2x2„-2Ci3-a:2)(2x^-2C'i4-a;2) 

a:^ (-3Ci3 /2 - 3C14/2 - ^34 - 6) + 2x^ 
(2a;^-2Ci3-a;2)(2a;^-2Ci4-a;2) 



(B15) 



-'da 



(1 



\3/2 



(i51idri3 



2C13CI4 +4C34 — 8C13C34 — 8C14C34 



+ 



+ 



(2a;^-2Ci3-a;2)(2a;^-2(:7i4-a;2) 

("^14 ~3C|4 + 6Ci3 + 8Ci4 + 2Ci3C34 + 2Ci4C34 + 8(734) 

(2x^-2C'i3-x2)(2x2^-2C'i4-a;2) 
x^(-3Ci3/2-2C34-9)+3a;^/2 



(2a;^-2Ci3-a:2)(2x^-2Ci4-x2) 



(B16) 



2Cl^ - 2(7|4 +4(713^34 +4(7i4(734 +ar^H(7i3 - 4(734) +4a:^ 
(2a;^-2(7i3-x2)2 



(B17) 



th23 

-"dn 



2(7i4(7|4 +4(7|4 — 8(7i3(734 — 8(7i4(734 



+ 



+ 



(2a;^-2Ci3-a;2)(2x^-2Ci4-x2) 

a^m(-<^13-3C'|4 + 10gl3 + Cl3(7l4 + 6Ci4 + 2C'i3C34 + 2(7i4C'34 + 8C34) 

(2x^-2(7i3-x2)(2a:^-2(7i4-x2) 

a:^(-^i3-2(7i4-2(734-10)+2a:^ 
{2x1 - 2(7i3 - ) (2a;^ - 2(7i4 - ar^ ) 



(B18) 



4(734 — 8(7i3(734 — 8(7i4(734+a:^(8(7i3 + 8(7i4 + 8(734) — Sxf^ 
(2a;^-2(7i3-a;2)(2a;^-2Ci4-a;2) 



(B19) 



I^n = (l-xlf jdn.d^l. 



-2Cl^ + 2(7i4(734 ((7i3 + 6(7i4 - 36^34 ) - 3x^ 
{2xl,-2CrA-xiy 



(B20) 



^^34 
-'dn 



-(1 



2 \3/2 



/ 



2(7^3 — 2(7i3(7i4 + 2(7i3(7i4 — 2(7i3C34+2(7i3(7i4(734 



(2a;^-2Ci4-x2)2 



+ 



+ 



(-^13 - + 2(7i3 +4^13^14 + 2(7i4) 

(2x^-2Ci4-x2)2 
xt(-2(7i3-2Ci4-2) + 24 



(2ar^-2Ci4-a;2)2 



(B21) 



rM44 
'dO 



fi 2 n3/2 / ,^ 2(7i^3 - 2(7|4 +4(7i3(734 + 4Ci4(734 H(7i4 - 4(734) +4a:^ 
[l-x^) /dl]idl]3 (2x^-2Cm-x?)2 



(B22) 
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